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ABSTRACT 

I ■ Aims. We present the results of a detailed analysis of 452 ground-based high-resolution high S/N spectroscopic measurements spread over 4.5 
' years for /3 Canis Majoris with the aim to determine the pulsational characteristics of this star, and to use them to derive seismic constraints on 
^v^J , the stellar parameters. 

£ — ■ Methods. We determine pulsation frequencies in the Si III 4553 A line with Fourier methods. We identify the m- value of the modes by taking 
into account the photometric identifications of the degrees C. To this end we use the moment method together with the amplitude and phase 
. variations across the line profile. The frequencies of the identified modes are used for a seismic interpretation of the structure of the star. 
' Results. We confirm the presence of the three pulsation frequencies already detected in previous photometric datasets: f\ = 3.9793 cd _1 
"^h (46.057 yuHz), f 2 = 3.9995 cd~' (46.291 yuHz) and f 3 = 4.1832 cd~' (48.417 yuHz). For the two modes with the highest amplitudes we 
Q ■ unambiguously identify (ti,tn\) = (2,2) and (ti, m%) = (0,0). We cannot conclude anything for the third mode identification, except that 
m 3 > 0. We also deduce an equatorial rotational velocity of 31 ±5 kms~' for the star. We show that the mode f t must be close to an avoided 
crossing. Constraints on the mass (13.5 ± 0.5M o ), age (12.4 ± 0.7 Myr) and core overshoot (0.20 ± 0.05 Hp) of /? CMa are obtained from 
seismic modelling using f\ and fi. 
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. 5~ . 1- Introduction reasonably evolved star or oscillates in modes different from 

the fundamental. 

Many breakthroughs have recently been achieved in the field 

of asteroseismology of (3 Cephei stars. The observation of a The variability of /? CMa has been known for on e century 

few pulsating modes led to constraints not only on global stel- and th e star has been exten s ively studied. W e refertojAlbrecht 

lar parameters but also on the core overshoot parameter and <ll908l>. iHenroteaul Jl918l). iMeveJ Jl934 and IStruvd Jl950l) 

on the non-rigid rotation of several f3 Cephei stars. In particu- for the first spec troscopic measurements of /3 CMa. Later 

lar modelling has been pe rformed for HP 129929 jAerts et all IShobbrookl Jl973l) found three pulsation frequencies from ex- 



2003a : iDuoret et alJl2004 and v Eri JPamvatnvkh et alJl2004t tensive photometric time s eries. The same t hree f requencies 
Ausseloos et al.ll2004 . Our aim is to add other ft Cephei stars were recently confirmed bvlShobbrook et all 0006) who anal- 



to the sample of those with asteroseismic constraints. ysed photometric measurements of a multisite campaign dedi- 

The B1II-III bright /3 Cephei star p Canis Majoris cated to the star. 
(HD 44743, HR2294, V mag = 1.97) is particularly interesting lAerts et all ill 994 coUected spectroscopic data in order to 
to study. Indeed, earlier photometric and spectroscopic data identify the modes of the known frequencies of ft CMa. In this 
revealed that this object exhibits multiperiodicity with rather papel . we present a similar ana i ysis but based on a much i^ger 
low frequencies in comparison with the frequencies of other numb erof s pectra and using the vers ion of the moment method 
P Cephei stars, which would indicate that p CMa is either a improV ed hv lBrimiet fc A~erts] J2003I) . We then construct stellar 
models which show oscillations in accordance with our unique 

Send offprint requests to: A. Mazumdar identification of the modes of /? Canis Majoris. 

* Based on spectroscopic data collected with the CORALIE echelle 
spectrograph attached to the 1.2m Swiss Euler telescope at La Silla, The paper is organised as follows. Section describes the 

Chile. results from our spectroscopic observations, including data re- 

** Postdoctoral Fellow of the Fund for Scientific Research, Flanders, duction, frequency analysis and mode identification. In Sect. [3] 
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Table 1. Observing logbook of our spectroscopic observations 
of /3 CMa. 



Number of JD 
observations 2450000 + 





Start 


End 


22 


1591 


1598 


26 


1654 


1660 


70 


1882 


1894 


62 


1940 


1953 


60 


2227 


2239 


27 


2569 


2582 


2 


2624 


2624 


91 


2983 


2994 


54 


3072 


3085 


38 


3271 


3283 



we present our seismic interpretation of (3 CMa. We end the 
paper with a discussion about our results in Sect.|4] 

2. Spectroscopic results 

2.1. Observations and data reduction 

Our spectroscopic data were obtained with the CORALIE 
echelle spectrograph attached to the 1.2m Leonard Euler tele- 
scope in La Silla (Chile). Because the beat period between the 
two known dominant frequencies is about 50 days we collected 
data on a long time span. Observations were collected during 
several runs spread over 4.5 years. The number of observations 
and the ranges of their Julian Dates are given in Table ^ in 
total, we gathered 452 spectra during 1692 days. 

An on-line reduction of the CORALIE spectra, using the 
INTER-TACOS software package, is av ailable. For a descrip - 
tion of this reduction process we refer to Barann e et alJ (11996). 
We did a more precise correction for the pixel-to-pixel sensi- 
tivity variations by using all available flatfields obtained dur- 
ing the night instead of using only one Hatfield, as is done by 
the on-line reduction procedure. Finally, all spectra were nor- 
malised to the continuum by a cubic spline function, and the 
heliocentric corrections were computed. For our study of the 
line-profile variability we used the Si III triplet around 4567 A. 
This triplet is very suitable to study ft Cephei stars since the 
lines are strong, dominated by tem perature broadening, an d not 
too much affected by blending (see lAerts & De Cal2003l) . 
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Fig. 1. From top to bottom: phase diagram of < v 1 > for ft = 
3.9793 c d -1 , for f 2 = 3.9995 c d" 1 after prewhitening with ft 
and for ft =4.1832cd~' after prewhitening with ft and ft. 



2.2. Frequency analysis 

We performed a frequency analysis on the first three yeloc- 
ity moments < v 1 >, < v 2 > and < v 3 > (see lAerts et alJI 19921 
for a definition of the moments of a line profil e) of the Si III 
4553 A line by means of the program Period04 ( Lenz & BregeJ 
20051). For some (i Cephei stars (se e Schriivers et al. 2004; 
Telling et alll997tlBriauet et al.l2005h a two-dimensional fre- 
quency analysis on the spectral lines led to additional frequen- 
cies compared to the one dimensional frequency search in inte- 
grated quantities such as moments. Consequently we also tried 



to find other frequencies for /3 CMa by means of this latter 
method. 

Despite the strong aliasing of our dataset, our frequency 
analysis reaffirmed the presence of the already known three 
pulsating frequencies of yS C Ma. ft = 3.9793 c d' 1 , f? = 3.9995 
c d" 1 and f 3 = 4.1832c d 1 I Sho bbrookll 9*731 IShobbrook et al J 
2006). Unfortunately no other frequencies could be discovered 
in our new spectroscopic data. Phase diagrams of < v 1 > for 
ft, ft and ft are shown in Fig.^ The frequencies and ampli- 
tudes that yielded to the best fit of < v 1 > are listed in Table|2] 



Mazumdar et al.: Asteroseismology of /3 Canis Majoris 



3 




40 60 



-60 -40 



-20 20 

V (km s" 1 ) 



40 60 




-60 -40 



-20 20 

V (km s" 1 ) 




-60 -40 -20 20 

V (km s" 1 ) 



Fig. 2. Amplitude and phase distributions for /i = 3.9793 cd , f 2 = 3.9995 cd~' and f 3 = 4.1832 cd 
line. The amplitudes are expressed in continuum units and the phases are expressed in n radians. 



for the Si III 4553 A 



Table 2. Frequencies and amplitudes of the first moment of 
the Si III 4553 A line together with their standard errors. The 
quoted errors for the frequency are intermediate between the 
overestimated value given by the frequency resolution, l/T ~ 
6 x 1CT 4 c d _I , and the underestimated value given by Period04 
for an ideal case free from aliasing and for white uncorrelated 
noise (6 x 1(T 6 ,6 x 1(T 6 and 2 x 1(T 5 cd" 1 for f u f 2 and f 3 
respectively). 







Frequency 


Amplitude 




(cd 




0"Hz) 


(kms" 1 ) 


f\ 


3.9793 ± 


0.0001 


46.057 ±0.001 


2.7 ±0.1 


h 


3.9995 ± 


0.0001 


46.291 ±0.001 


2.6 ±0.1 


h 


4.1832 ± 


0.0001 


48.417 ±0.001 


0.7 ±0.1 



These three frequencies reduce the standard deviation of the 
first moment by 72%. Fig. |2] shows the variations across the 
Si III 4553 A line for f, , h and / 3 (see iTelting & Schrijversl 
ll997tlSchriivers et all 199^ for a definition). 

We also performed a frequency analysis on the equivalent 
width but we could not find any significant frequency. This 
is not uncommon since significant EW variations have been 
found in o nly a few Cephei star s with photometric variations 
up to now JDe Ridder et alJ20o3l . 

2.3. Mode identification 

Our methodolo gy to identify the mo des of p CMa is similar to 
the one used in lBriauet et al.l J2005h . which led to a successful 
mode identification for the p Cephei star 9 Ophiuchi. We refer 
to that paper for a detailed explanation on our chosen process. 



2.3.1. Adopted ^-identifications 

We make use of spectroscopy to identify the values of the az- 
imuthal number, m, which are not accessible from photometry. 
In order to limit the number of parameters in our spectroscopic 
mode identification we adopt the degree, (, as obtained from 
photo metric mode identification. Recently, Shob brook et alJ 
(2006) found t\ - 2 for the first mode, f\, and l 2 = for the 
second mode, f 2 . However, the degree of the third mode could 
not be determined. 

With our dataset we corroborate the mode wi th f 2 
to be radial as foll ows. [Telting & Schriiversl Jl997l) and 
ISchriivers etal] i ll 9971) showed that, when there is a minimum 
(almost zero) in the amplitude and a corresponding phase shift 
of n near the centre of the line profile, one can conclude to be 
dealing with a radial or a dipole mode. As evident from Fig. [2J 
this is clearly the case for f 2 . 

2.3.2. Moment method 

With the adopted ^-values we then determine the m-values 
by means of the moment method. The new implementation 
of this techniqu e was o ptimised for multiperiodic signals by 
Briq uet & Aertsl J2003). Its improvements and the huge in- 
crease in dataset explain why w e obtain a mode identification 
different from Aert s et alJ ([1994) who used an older version of 
this method and whose data did not cover the beat periods of 
P CMa. 

The theoretical moment values to be compared to observed 
moment values are computed by fixing the following param- 
eters. A linea r limb-darkening co efficient u of 0.292 is taken 
(see e.g., Wad e & Rucinskil ll985). The ratio of the amplitude 
of the horizontal to the vertical motion, denoted by K, is given 
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Fig. 3. The theoretical amplitude distribution for f\ computed 
from line profile time series generated for the best parameter 
combinations derived with the moment method. 



by K — GM/a> 2 R 3 , where M is the mass, R the radius and to the 
angular pulsation frequency. With M = 13.5M and/? = 7 &R & 
foe Caa2002llHevnderickx et alJl99i) we obtain Ki = 0.133, 
K2 =0.132 and ^3 = 0.121. We varied the free parameters in 
the following way: the projected rotation velocity, v sin i, from 
1 to 35 km s with a step 1 km s , the inclination angle of the 
star, i, from 1° to 90° with a step 1°, and the line-profile width 
due to thermal broadening, cr, from 1 to 20 kms -1 with a step 
1 kms . 

The mode identification by means of the moment method 
gives a preference to nt\ = 2 but we cannot rule firmly out nt\ 
I . For the mode with frequency we cann ot conclude any- 
thing , as is the case for the photometric data JShobbrook et alJ 

2.3.3. The amplitude and phase variations across the 
line profile 

Because the discriminant values of the best moment solutions 
are very similar we need an additional check in order to safely 
conclude that (i\,m\) = (2,2). Our method is to visualise the 
behaviour of observed amplitude and phase variations across 



the line profile compared to theoretically computed ones for 
the best parameter sets given by the moment method. 

The observed amplitude and phase variations across the 
Si III 4553 A line are shown in Fig.[5]for f x , f 2 and fa . The theo- 
retical distributions were computed from line profile time series 
generated by means of Townsend's codes llTownsendlll997h . 
called BRUCE and KYLIE. The line-profile variations as well 
as the amplitude and phase variations were computed by con- 
sidering the three modes together. 

The amplitude distributions for f\ computed for the best pa- 
rameter combinations for (€\,m{) = (2, 2) and (t\,m{) = (2, 1) 
given by the moment method are shown in Fig. [3] By com- 
paring the observed (see Fig. |3 to the best theoretical ampli- 
tude distributions (see Fig. [3 for the mode with frequency fi, 
it becomes clear that ((\,m\) = (2,2). Indeed the behaviour 
of the theoretical amplitude distribution differs from the ob- 
served one for all best solutions with (€\,m{) = (2, 1) since the 
"triple-humped" character of those solutions are absent in the 
observed amplitude distribution (see top left panel of Fig. |2j. 
All best parameter combinations with (€\,m{) - (2,2) mimic 
very well the observations. For the third mode we unfortunately 
could not discriminate between the different solutions. All we 
could derive from its phase behaviour in Fig.|2]is that 1113 > 0. 

2.4. Derivation of the stellar equatorial rotational 
velocity 

Each solution given by the moment method indirectly gives a 
value for the equatorial rotational velocity, since the inclina- 
tion i and the projected equatorial velocity vq = v eq sin i are 
estimated. We made a histogram (see Fig. |4j for v eq by con- 
sidering only solutions with (i\,m\) = (2,2) and by giving 
each equatorial rotational velocity v^,* its appropriate weight 
wit = Eo/Efc, where So is the discriminant value for the best so- 
lution. By calculating a weighted mean and standard deviation 
of the data we got v eq = 31 ±5 kms"' for/3 Canis Majoris. We 
also constructed a histogram for the inclination angle in a sim- 
ilar way. We found a flat distribution in the range [0°,90°] so 
that we could not restrict the value for 2. The moment method 
could consequently allow us to limit the range for the couple 
(v sin i, i) but not v sin i or i separately. 

3. Seismic interpretation 

We have performed a thorough seismic analysis of the observed 
frequencies by comparing them with frequencies obtained from 
theoretical stellar models. Apart from the two identified fre- 
quencies, f\ and fi, we have also taken into account the posi- 
tion of fi CMa on the HR diagram in our modelling. The third 
unidentified frequency, fa, is not used in the modelling process. 

3.1. Global parameters 

The effective temperature of /3 CMa has been quoted in the lit- 
erature in a range of val ues aroun d 24000 K. The extreme val- 
ues a re log = 4.37 (Hevnderickx et al. 1994) and log r e ff — 
4.45 dTian et alJ2003h . Recentlv lMorel et alJd2006l) made a de- 
tailed NLTE analysis of /3 CMa to find log r eff = 4.38 + 0.02. 
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Fig. 4. Histogram for the equatorial rotational velocity of the 
star. 



We shall adopt a conservative range of log T e ff = 4.40 + 0.04 
for the present work. 

The luminosity can be determine d from t he Hipparcos par- 
allax: n - 6.53 ± 0.66 mas (Perr vman et al.ll 19971) . Assuming 
a bolometric correction of -2.29 mag, this translates to a lumi- 
nosity range of log L/L Q = 4.41 + 0.16. However, the value 
of the bolometric correction might be a major source of er- 
ror in this calibration. We also fin d a range of values in the 
literature, from logL/ L = 4.40 dStankov & Handler! Eo05) 
to log L/L Q = 4.79 dTianetalJ|2003l) . We adopt the range 
log L/L Q = 4.45 + 0.20, which covers most of these quoted 
values. 

The metailicity of j3 CMa has been found 

to be [M/H] = 0.04 + 0.10 in recent studies 

(N iemczura & Daszvnska-Daszkiewiczl l2005h . We have 
adopted this range in our models. Assuming a solar metailicity 
of Z = 0.018, this translates into a range of metailicity for 
P CMa: 0.016 < Z < 0.023. For a given metailicity, we have 
also varied the hydrogen abundance slightly while exploring 
the parameter space for the models. 

3.2. Rotational splitting 

Since the rotational velocity of [3 CMa has been found to be 
moderately low, we do not expect any strong coupling between 
rotation and oscillation frequencies. However, one cannot be 
certain if the interior rotation velocity is higher or not. Indeed, 
evidence of differential r otation of the core has been found in 
two o ther [3 Cephei stars (tPupret et al.l2004llPamvatnvkh et all 
120041) . But in the absence of observations of a rotational mul- 
tiplet, or a definite estimate of the rotation velocity, we do not 
have enough constraints to check if the internal rotation is in- 
deed uniform or not. As a first approximation, therefore, we 
have assumed rigid rotation in the interior of the star. 

We constructed non-rotating stellar models and accounted 
for the effect of rotation by adding the rotational splitting 
6v s - m/3v eq /R to the theoretical frequency, where j3 represents 
the Ledoux constant and depends on the rotational kernels. The 



value of f3 is calculated for each mode from the eigenfunction. 
The rotational splitting is found to lie between 0.9 /iHz and 
1.3 //Hz, depending on the model, the uncertainty being pri- 
marily due to the rotational velocity, v eq . Since 6v s depends on 
/3, instead of using a uniform value for 6v s , we have calculated it 
for each particular mode of every model that we want to com- 
pare with the observations. Thus, we calculate the theoretical 
value of the rotationally split component (I = 2, m = 2) of the 
model frequency: v mo d e i(n, t — 2, m — 2) = v mot j e i(n, € = 2, m = 
0) + 6v s and match it with the observed frequency f\ . 

3.3. Stellar models 

We constructed a grid of stellar models between masses 
of 12Mp, and 16M Q using the CESAM evolutionary code 
dMorell Il997l). These models use d the OPAL equation of 
state ( Rogers & Navfonov 2002) and OPAL opacity ta- 
bles (Iglesias & Rogers! Il996l). complemented by the low - 
temperature opacity tables of [Alexander & Fergusonl dl994l) . 
Con vection was describ ed by the standard mixing length the- 
ory dHenvev et alJll965l) and nuclear reaction rates wer e ob- 
tained from the NACRE compilation dAngulo et all 19991) . The 
frequencies of oscillation were computed under the adiabatic 
appr oximation, using the Aarhus puls ation package, ADIPLS 
dChristensen-Dalsgaard & BerthomieJl99ll). 

Our models adequately span the range of stellar parame- 
ters such as mass (12 < M/M Q < 16), chemical composition 
(-0.06 < [M/H] < 0.14) and core overshoot (0 < d ov < 0.3, 
where d m represents the extent of overshoot in terms of the 
local pressure scale height). For each combination of these pa- 
rameters we constructed models spanning the main sequence 
phase of evolution. The frequencies of models lying within or 
close to the error box on the HR diagram were then compared 
to the observed frequencies of CMa. The mixing length of 
convection was not varied since it does not play a significant 
role in the models, the outer convective envelope being ex- 
tremely thin in such stars. Diffusion and radiative levitation of 
elements were not incorporated in the models. 

3.4. Matching the frequencies 

Although the (£, m) values of two of the modes were identi- 
fied from their photometric and spectroscopic characteristics, 
their radial orders cannot be determined from observations. 
Therefore, we have to account for various possibilities for the 
radial order in the seismic modelling. First we deal with the ra- 
dial mode, f%, by comparing it with the radial frequencies of the 
stellar models. It becomes evident that the radial mode must be 
either the fundamental mode (( — 0, n — 1 ) or the first overtone 
({ = 0, n = 2) for the star to be at its prescribed position on the 
HR diagram. If fz were to be any higher order radial frequency, 
then the radius of the star (hence the luminosity and T e ff) would 
become completely inconsistent with the estimates of log T e ff 
and log L/L G for /3 CMa obtained independently. Indeed, ei- 
ther of the fundamental mode or the radial overt one has b een 
observed in other 13 Cephei stars as well (e.g ., iDupret et alJ 
120041 lAusseloos et ail [2004: Ae rts et alJExxl . We explored 
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both possibilities by comparing the observed frequency, fa with 
the theoretical radial mode frequencies of our models. 

3.5. Radial mode - fundamental or overtone? 

As a first step, we matched fa to the radial fundamental mode 
of our models. The frequency fa was then compared to differ- 
ent orders of Z — 2 g- and /?-modes (corrected for rotational 
splitting, as explained in Sect. 13.21 It turns out that for the fre- 
quency range of interest, this mode being the first g-mode (g\: 
Z — 2,n — -1) is the only possible solution. The so-called f- 
mode (Z = 2, n = 0) is too high and the second g-mode (g2- 
Z = 2, n = -2) is too low in value compared to fa - 6v s . 

However, although it is easy to find a number of mod- 
els with different stellar parameters whose radial fundamental 
mode matches with fa, none of these models match the fre- 
quency fa . The central problem lies in the fact that the g\ mode 
frequency of any given model is too close to the radial funda- 
mental mode, compared to the distance between them as indi- 
cated from the observations of fa and fa. The only other nearby 
Z = 2 mode, the g2 mode, is far too distant from the radial fun- 
damental. In other words, the "separation" between the radial 
fundamental mode and the 1-2 modes in the models is ei- 
ther much smaller (for the g\ mode) or much larger (for the g2 
mode) than the observed value, fa - (fa - 6v s ). What we require 
to match both frequencies is an intermediate value (~ 1.2 yuHz) 
of the separation between the radial mode and the Z = 2 mode. 

The clue to our problem lies in recognising the fact that 
such a small (but not too small) difference between the £ = 
and 1-2 modes can only occur in case of an avoided crossing. 
As long as the modes follow the regular smooth separation pat- 
terns, we will never find a solution where the Z — and 1-2 
modes have the appropriate separation as required by the ob- 
servations. This is illustrated in Fig. [5] In this figure, we show 
how the Z — and i = 2 modes vary with evolution of a given 
star. For the purpose of illustration, we have assumed an aver- 
age value of the rotational splitting (with an error margin) and 
corrected the observed frequency, fa, to compare it with the 
theoretical frequency. As the star evolves, the density gradi- 
ent at the edge of the shrinking convective core increases. This 
has the effect of rapidly increasing the frequencies of the g- 
modes, which causes successive "bumping" of the modes. For 
example, as shown in Fig.|5j the g2 mode of the 1 = 2 degree 
first increases (for X c > 0.16), until it bumps into the g\ mode 
(around X c 0.16), creating an avoided crossing between these 
two modes at that age. With further evolution, the gi mode in- 
creases in value until it creates another avoided crossing with 
the /-mode (around X c =s 0.12). Progressively, higher order 
modes are bumped as the star evolves. An excellent discussio n 
of such mode bumping is provided bv lAizenman et al] Jl977l) . 

We find that as long as we are constrained by the position of 
/? CMa on the HR diagram, the required separation between the 
Z = and Z = 2 modes (~ 1.2 yuHz around a frequency value of 
45 //Hz) cannot be obtained if the radial mode is the fundamen- 
tal mode. Even for an avoided crossing, the appropriate sepa- 
ration between the radial fundamental and the Z = 2,n = — 2 
bumped up mode occurs at a much lower frequency than ob- 
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Fig. 5. The variation of the frequency of Z = and 1 = 2 
modes of a 13M stellar model with evolution (characterised 
by the central hydrogen abundance, X c ) is shown. The dotted 
lines indicate the Z = modes, while the dashed lines show the 
Z = 2 modes. Each mode is labelled by its degree and radial 
order (Z, n). The horizontal solid line and the shaded strip show 
the observed radial mode (fa) and the splitting-corrected Z = 2 
mode (fa —6v s ) respectively. The vertical dot-dashed line shows 
a possible solution where the theoretical frequencies match the 
respective observed frequencies. 



served. Fig. |5J which represents a typical model at the appro- 
priate position on the HR diagram, confirms this. However, a 
solution is indeed plausible if we consider the observed radial 
mode to be the first overtone instead. In that case, the avoided 
crossing between the Z = 2 g\ and /-modes bumps the g\ close 
to the Z = 0, n = 2 mode. We now, therefore, explore the possi- 
bility of the radial mode being the first overtone mode. 

By considering the radial mode as the first overtone (Z = 
0, n = 2), we do indeed find a number of models with dif- 
ferent stellar parameters that match the observed frequencies. 
Specifically, fa is matched with the Z = 0, n = 2 mode, and fa 
with the Z = 2,n = — 1 mode, after correcting for the rotational 
splitting. We have also checked the excitation rates of th ese two 
frequ encies by the nonadiabatic oscillation code, MAD ( D upretl 
2001) and found both of them to be excited for all the models 
which fit the data. 

A similar solution with the radial mode being the second 
overtone (Z = 0, n = 3) and the Z = 2 mode being the bumped- 
up /-mode is again ruled out because the relevant avoided 
crossing occurs at a higher frequency than observed. 

The common feature of all the possible solutions is the fact 
that the Z = 2,n = — 1 mode is an avoided crossing with the 
/-mode. This helps to constrain the stellar parameters a great 
deal, especially the age (or X c ) of the model. The requirement 
that the star must be at the precise age for a specific avoided 
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Fig. 6. The variation of the model frequencies with three stellar parameters - M, d ov and Zq - are shown. In each panel, two of the 
parameters are fixed to clearly bring out the effect of the remaining parameter on the frequencies. The I — 0, n — 2 frequencies are 
plotted as abcissae, while the I = 2, n = — 1 frequency, corrected for rotational splitting, are plotted as ordinates. The dashed lines 
indicate the two observed frequencies. The point where the two dashed lines intersect indicate the perfect match of the desired 
model. Each shaded band represents an evolutionary track. The width of the band along the y-axis indicates the uncertainty due 
to the rotational splitting. A full-colour version of this figure is available in the online edition. 



crossing to occur, is indeed a very strict one. Nevertheless, 
since we have only two frequencies to constrain our model 
with, we do find multiple combinations of the stellar param- 
eters where such a solution exists. We have investigated the 
limits of the stellar parameters that can be constrained with the 
observed frequencies (see next section). Although the limits on 
effective temperature and luminosity were not explicitly used 
in determining the best models, all of them do lie well within 
the errorbox on the HR diagram adopted in Sect. 13. II 

3.6. Limits on stellar parameters 

We have, in principle, four major stellar parameters to tune our 
models with - M, d ov , Xq and Zq. For each appropriate combi- 
nation of these parameters, the age (in terms of X c ) is automati- 
cally chosen by the closest match of the frequencies. The initial 
hydrogen abundance, Xo, must, however, be linked somewhat 
to the metallicity, Zq; we have varied Xo only within permissi- 
ble limits for a particular choice of Zq so that [M/H] remains 
within the bounds quoted in Sect. 13. II In Fig.|6] we have shown 
how the relevant frequencies of the models, v(i = 0, n = 2) 
and v(( = 2,n = -1), vary as each of the stellar parameters 
are changed. In each panel, only one parameter is changed at a 
time, keeping the others constant. Each band is, in fact, an evo- 
lutionary track - the sense of evolution being from the right to 
the left (decreasing radial frequency). The width of the shaded 
bands along the y-axis reflects the uncertainty due to the rota- 
tional splitting correction. The tracks are truncated beyond the 
ages where the models lie outside the error box on the HR di- 
agram. Our best models lie near the centre of each plot, where 
both the model frequencies match the observed ones. 

In Fig.[7]we vary multiple parameters simultaneously to ob- 
tain a good match of the model frequencies with the observed 
ones. This illustrates that several solutions are possible in the 
multi-dimensional parameter space. Actually, we have shown 
only the models with the extreme limits of the stellar param- 



eters for which a solution can still be found. These are only 
indicative of the trend of the frequencies of the models, and 
several other solutions are possible when the parameters (espe- 
cially overshoot and metallicity) are varied within their bounds. 
Table[3]lists the physical parameters for these selected models. 

We find that no solutions are possible in the absence of core 
overshoot. Even with a small amount of overshoot (d m = 0.1), 
we need to have a higher mass to obtain a solution. The higher 
stellar mass helps to increase the mass of the convective core to 
offset the low overshoot. Most of our best fit models use d ov = 
0.2. Higher overshoot (d m = 0.3) models can reproduce the 
frequencies only for low metallicity and proportionally lower 
helium content, again suggesting a trade-off between overshoot 
and helium content to maintain a balance of the core size. This 
indicates that the size of the convective core plays a crucial role 
in determining the frequencies. Indeed, all our best fit models, 
including those with low and high overshoot have fractional 
convective core mass of 24-28%. 

The mass of our best fit models are mostly limited to 
13.5 + 0.5M Q , for d m = 0.2. The mass could be higher (up to 
14.5M Q ) if either the metallicity is low ([M/H] = -0.06) or the 
overshoot is low (d ov = 0.1). The central hydrogen abundance 
of the best models are limited to the range 0.128 > X c > 0.123, 
for overshoot values of d m = 0.2. As expected, the low over- 
shoot models have younger ages and lower X c . The situation is 
opposite for high overshoot models. 

4. Conclusions 

/3 Canis Majoris is one of the f3 Cephei stars whose variabil- 
ity has been observed and analysed for one century. It was 
discovered that this star pulsates with three frequencies rather 
low in comparison with other known stars of its type. For this 
reason /? CMa is an important target for asteroseismology pur- 
poses. However, so far, no definite mode identification had been 
achieved for this star so that no modelling could be attempted. 
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Table 3. Physical parameters of representative stellar models that match the observed frequencies. 
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Fig. 7. Similar to Fig. [6] except that multiple parameters are 
varied at a time to obtain a good match between the theoret- 
ical frequencies and the observed ones. Selected models with 
limiting values of overshoot and metallicity are shown only. A 
full-colour version of this figure is available in the online edi- 
tion. 



Our aim was to increase the number of known pulsating fre- 
quencies and mostly to provide a unique identification of the 
modes of ft CMa. 

Our study was based on 452 ground-based high-resolution 
high S/N spectroscopic measurements spread over 4.5 years. 
We used the Si III 4553 A line to derive the pulsation char- 
acteristics of p CMa. Our dataset unfortunately suffers from 
strong aliasing but the three established frequencies of the star 
were confirmed in the first three velocity moments of the line 
and in the spectra themselves. They are fi = 3.9793 cd 1 
(fx = 46.057 /Mz), f 2 = 3.9995 cd" 1 (/ 2 = 46.291 /Mz) 
and f 3 = 4.1832 cd" 1 (/ 3 = 48.417 /Mz). Unfortunately no 
new frequencies were discovered neither in our spectroscopic 
observations nor in the recent mul tisite photometric measure- 
ments led bv lShobbrook et all J2006I) . 

The important result of the combination of both intensive 
campaigns is an identification of the two main modes of f3 CMa, 
which is a strong constraint for further asteros eismic m odelling 
of the star. The photometric identification by Shobb rook et alJ 
(2006) yielded t\—2 and i 2 — 0. Our spectroscopic data could 
corroborate that the mode with f 2 is radial. We adopted the 



photometric identification of £\ and spectroscopic techniques 
allowed us to derive the m- value of the main mode. The ap- 
plication of the moment method gave a preference to m\ = 2. 
Because moment solutions could not definitely exclude ni\ — 1 
we made use of the behaviour of the amplitude distributions 
across the line profile for the best parameter sets given by the 
moment method. In this way we could conclude without any 
doubt that (l\,m\) = (2,2). For the third mode nothing could 
be concluded. In addition we derived a stellar equatorial rota- 
tional velocity of 3 1 + 5km s . 

The definite identification of two of the observed modes 
and a much improved estimate of the rotation velocity of 
/3 CMa allowed us to attempt the first seismic modelling of this 
star. Although it is not realistic to hope for a unique model to 
fit just two frequencies, we have thoroughly explored the stellar 
parameter space to derive reasonable constraints for the mass, 
age, and core overshoot. The most significant aspect of the seis- 
mic analysis is the fact that we could assert that the non-radial 
mode, /i, is close to an avoided crossing. This implies a very 
strong constraint on the stellar parameters, especially the age of 
the star. At the same time, it rules out the possibility of the ra- 
dial mode, f 2 , being the fundamental mode. This makes /3 CMa 
one more (3 Cephei star kno wn to have a do minant radial over- 
tone mode of pulsation (cf. lAerts et all20 06). 

Our best fit models indicate that /3 CMa has a mass of 
13.5 + O.5M , an age of 12.4 + 0.7 Myr (X c = 0.126 + 0.003) 
and core overshoot of d m = 0.20 ± 0.05. No satisfactory model 
can be found if core overshoot is absent. A small overshoot 
parameter is possible only for a higher mass along with high 
metallicity (and proportionally higher helium content). On the 
other hand, higher core overshoot is required if the star is metal- 
poor ([M/H] < -0.05). However. lMorel etail fcOOfil) found the 
composition of CMa not much different from that of the Sun, 
making such possibilities unlikely. Therefore, it is safe to con- 
clude that the models with d ov = 0.20 + 0.05 are the most likely 
ones. If the chemical composition of /3 CMa could be known 
to higher accuracy independently, one would be able to con- 
strain the other parameters even more. All the solutions turn 
out to have effective temperatures close to the cooler edge of 
the adopted errorbox on the HR diag ram. This is cons istent 
with the recent estimates of T e g (e.g., M orel et al.l l2006). We 
note that we had deliberately chosen a very conservative error- 
box for effective temperature and luminosity; a stricter limit on 
these parameters would rule out some of our possible models. 

In retrospect, one can also try to identify the mode of os- 
cillation for the third frequency, fa, by comparing it with the 
theoretical model frequencies. In this comparison, we allowed 



Mazumdar et al.: Asteroseismology of /3 Canis Majoris 



9 



for different rotational splitting values, m^, for each non-radial 
mode with the restriction > 0, as we found in Sect. [2 This 
leads us to only one possibility for f$, for all the models which 
match f\ and fa, Co, = 3,ni3 = 2, no, = -1. We hope that this 
identification of fj can be checked through future observations. 
However, we cannot place further constraints on the models at 
this stage using this frequency; all the models with stellar pa- 
rameters in the range restricted by the first two frequencies also 
match the third frequency with the identification given above 
within the uncertainty associated with the rotational velocity. 
One would need a more precise estimate of the rotational ve- 
locity to distinguish between these models. 

The rotation period may be calculated from our estimate 
of the equatorial velocity (Sect. \2Al and the radius of our best 
models which indeed lie in a narrow range (Table [3j. We esti- 
mate the rotation period to be 18.6 + 3.3 days, which indicates 
that /3 CMa is indeed a slow rotator; therefore, our assumption 
in neglecting higher order terms of the rotation velocity while 
calculating the frequency splitting stands justified. 

Despite the knowledge of only two frequencies for this star, 
the occurrence of the avoided crossing goes a long way towards 
constraining most of the stellar parameters. While one cannot 
expect to be so lucky for every star, we have shown that the 
identification of an avoided crossing might help us to extract a 
lot more information about the star than any normal mode. 
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